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Abstract 

The aim of this paper is to investigate the distribution of a continuous homopoly- 
mer in the presence of an attractive finitely supported potential. The most intricate 
behavior can be observed if we simultaneously vary two parameters: the temper- 
ature, which approaches the critical value, and the length of the polymer, which 
tends to infinity. As the main result, we identify the distributions that appear in 
the limit (after a diffusive scaling of the original polymer measures) and depend on 
the relation between the two parameters. 

1 Introduction 

We consider the following model of long homogeneous polymer chains in an attractive po- 
tential field. Let the space C([0,T],M d ) be equipped with the Wiener measure Pot and 
suppose that we have a smooth, nonnegative, not identically equal to zero, compactly sup- 
ported potential v G C x> (R d ) and a coupling constant (3 > (inverse temperature), which 
regulates the strength of the attraction. The elements of the space are interpreted as 
realizations of a continuous homopolymer on [0,T] and are distributed according to the 
Gibbs measure Pp,T with 

dPp,T, , ePtfvHWt 



u) = , weC([0,n: 



where 

Zp, T = Eo tT e?fi v( - u ® )dt 

is the partition function. 

Thus the polymer measure we consider is of a "mean field" type, where the polymer 
chain interacts with the external attractive potential (as in [6]). While the potential is 
assumed to be constant in time in our paper, many interesting results have been ob- 
tained for disordered media, i.e., time-dependent random potentials (existence of phase 
transitions, dependence of the growth rate of the partition function on the temperature, 
etc., see [2J, [5], for example). We suspect, however, that the disordered models are not 
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amenable to such a detailed analysis of the phase transition phenomena as found in the 
current paper. 

It is established - using the Feynman-Kac formula - that under the measures Pp,r, the 
processes {u(t),0 <t < T} are time-inhomogeneous and Markovian and that their tran- 
sition densities can be expressed in terms of the fundamental solution pp of the parabolic 
equation 

j^ = Hpu, where Hp = -A + (3v : L 2 (R d ) — )■ L 2 (R d ). (1) 
More precisely, the finite-dimensional distributions are 

P PiT (u(ti) e A u u{t n ) G A n ) = (2) 
= — — / ... / / p /3 {t 1 ,0,x 1 )...p p {T - t n ,x n ,y)dydx n ...dx 1 

for < ti < ... < tn < T and A u A n G B{R 3 ). 
It is worth noting that 



Z/3,t = / pp{t,Q,y)dy. 



It was shown that in d > 3, at a ceratin (critical) value of the coupling constant, there 
occurs a phase transition occurs between a densely packed globular state and an extended 
phase of the polymer. Namely, for (3 > /3 cr , a typical polymer realization is at a distance 
of order one from the origin as T — > oo. On the other hand, for /3 < /3 cr , the realizations 
need to be scaled by the factor \JT in the spatial variables in order to get a non-trivial 
limit (as in Theorem II .11 below). The critical value of the coupling constant is 

/3 cr = sup{/3 > 0| sup cr(Hp) = 0}, 

where cr(Hp) is the spectrum of the operator H^. 

Here's the precise statement of the result proved in [3] that is relevant to this paper. 

Theorem 1.1. For a > 0, let f a : C([0,T],M 3 ) -»■ C([0, aT], M 3 ) be the mapping defined 
by (f a u)(t) = \fauj{t/a), and let f a P be the push-forward of a measure P by this mapping. 
There is a measure Qo on C([0, 1], R. 3 ) that corresponds to a certain time-inhomogeneous 
Markov process starting at the origin such that 

/T-iPjSor.T =>- Qo as T -> oo. 

Remark 1.2. We use the subscript in the notation for the limiting measure since we 
will introduce a whole class of measures, and Qo will be just a particular member of this 
class. 

Remark 1.3. For (3 < (3 CT , the measures fr~ l Pp,T converge to the Wiener measure on 
C([0, 1], M 3 ). Note, however, that this convergence and the convergence in Theorem 11.11 
take place for fixed values of (3 that don't scale with T. 
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In this paper we'll allow /3 = /3(T) to change in such a way that (/3(T) — /3 CI )yT — > x- 
We'll show that there are measures Q 7 and a linear mapping j(x) = c-X such that 

f T -iPp(r),T => Q 7 ( x ) as T ~> 00 . 

Before formulating this as a theorem, let us describe the limiting measures Q 7 . They 
were introduced in [4] as the polymer measures on C([0, 1], 1R 3 ) corresponding to zero- 
range attracting potentials (i.e., the potentials that are, roughly speaking, concentrated 
at the origin). More precisely, it was shown in |4] that, for each 7 G I, the polymer 
measures corresponding to the potentials 

converge as e \. 0. (Here xb is the indicator function of the unit ball centered at the 
origin.) The limit will be denoted by Q 7 . It is worth noting that all the measures Q 1 are 
spherically symmetric and thus the dependence of the polymer measures Ppcr),T on the 
'shape' of the potential disappears in the limit. 

Another, less intuitive, but more convenient, way to define the measures Q 7 is through 
the self-adjoint extensions of the Laplacian acting on IR 3 \ 0. Namely, it was shown in 
[lj that there is a one-parameter family {£ 7 ,7 G M} of self-adjoint operators acting on 
L 2 (R 3 ) such that C y f = Af whenever / G C °°(M 3 \ 0). The kernel of exp(t£ 7 ), t > 0, is 
given by 

e -\x-y\ 2 /2t l r e -V2X(\x\ + \y\)+Xt 

p (t, x, y) = + -— / — = r-~-^d\, (4) 

(2-Ktyi 2 Airh J T{a) (V2A - i)\x\\y\ 

where x, y 7^ 0, T(a) = {z G C|Re2; = a}, and a > 7 2 /2. Thus p (t,x,y) can be 
interpreted as the formal fundamental solution of the parabolic equation 

du 

- = C lU . 

By analogy with (j2J), we can define the measures P 7j t5 whose finite-dimensional distribu- 
tions are given by the formula 



Ej iT (w(*i) G A h ...,u{t n ) G A n ) = 

p 7 (t x , x, Xi)...p 7 (T - t n , x n , y)dydx n ...dx 1 



7>- 

1 

Z 1>T {x) J A! JA n 

for < ti < ... < t n < T and A u A n G S(R 3 ), where 

Z iA x )=l Pi{t,x,y)dy = l + — — = TT-r dX - 

Note the dependence of the measures on the initial point x G M 3 . In fact, neither 
p 7 (ti,x, y) nor Z 1 ^{x) are defined for a; = 0, but we can make sense out of P 7)T by 



taking the limit of P T as x — > 0. If we put Q 7 = P 1; we'll obtain the same family of 
measures corresponding to zero-range potentials. 

It is also worth mentioning that the assumption d = 3 is important. Indeed, there is 
no phase transition for d = 1,2, while for dimensions higher than 3 there are no nontrivial 
self-adjoint extensions of the Laplacian. 



2 The main result 

It was shown in Lemmas 5.4 and 5.5 of [3] that the solution space of the problem 

1 dib 

-Aijj + (3 cr v(x)ip = 0, Mx) = 0(\x\ --(x) = 0(\x\- 2 ) as Id oo 

2 or 

is one-dimensional, and ip can be chosen to be positive. 
The main result of our paper is the following. 

Theorem 2.1. If (/3(T) - f3 CT )VT ->• x, then 

/t-iP/3(t),t Q 7 ( X ) as T ^ oo, 

where -y(x) — C X w ^h 

(/ R 3 v(x)ij(x)dx) 2 



c= V2/(/3 c 2 rTl ) and 7 i 



A/27T J R3 f (x)?/' 2 (a;)(ix 



The proof will rely on the asymptotic formulas for the fundamental solution of ([T]) and 
the partition function, which can be found in the right-hand side of (T5]). We formulate 
these asymptotic formulas in several propositions below. 

Let us introduce the space 



C cxp (R") = / G C( 



c cxp := sup(|/(x)|e |:c|2 ) < oo 



For / G C exp (]R 3 ), the solution of the parabolic problem 



du 
~dt 

is given by the inverse Laplace transform 



— = Hp(T)U, u(0,x) = f(x) (5) 



u pm (t,x) = ~ [ e xt (Rp {T) (\)f)(x)d\, (6) 

Zm Jr{Xo(fi(T))+6/T) 

where A (/3) = supcr(if / 3), 5 > 0, and -Rg(A) = (ff^ — A/) -1 is the resolvent of Hp. 
Proposition 2.2. For e > 0, we have 
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1. (The solution of the Cauchy problem) 

1 a(f) f e# e'^^T 

with 7 = A/2x/(7i/?cr); ^ > 0, and 

where the error term satisfies 

sup _|g / (T,t,x)| = ||/|| Ccxp O(T- 3 / 2 ), (8) 

eT<t<T,eVT<\x\<e- 1 VT 

2. (The fundamental solution) 

1 Kip(y) f e^ e~^^^r 
p m (t, y, x) = —— M rfA + <KT, t, », x), (9) 

where 

_ sup T|g(T,t,w,x)| ->■ 0. 

We'll also need the asymptotics of the fundamental solution when y is of order y/T. 
Proposition 2.3. For e > 0, 

At . att M + H 

11/" e" 5 " vr 

Pp(T)(t,y,x) = p (t,y,x) + -= — / — = — — d\ + q(T,t,y,x), (10) 

VT47T i J T (£+s) (V2A - 7)|x||y| 

where po is the fundamental solution of the heat equation and 

sup T 3/2 q{T,t,y,x) ->■ 0. (11) 

eT<t<T,ev / T<|y| 1 M<<:- 1 V / T 

We will need one more result concerning the behavior of the partition sum with respect 
to the measure Pg T (the Wiener measure P ,t translated by the vector x). 

Proposition 2.4. If Zp jt (x) = ~E^ T e^ v ^ t )) dt } then for every e > 0, we have 
w/iere 

lim sup |g z (T, t, x)| = 0. 



3 Proofs 



Note that the expression ([6]) for the solution of (JSJ) contains the resolvent of the operator 
Hp inside the integral. It will be seen that the main contribution to the integral comes 
from the values of A that are close to zero. Therefore, it is important to know the 
asymptotics of the resolvent as A — > 0. First, let us make several observations concerning 
the spectrum of H p (a more detailed discussion of the spectral properties of Hp together 
with the proofs of the following four lemmas can be found in [3]). 
It is well known that for some < iV < oo 

a(^) = (-oo,0]U{A J }f =0 , 

where the eigenvalues are enumerated in decreasing order. 

Lemma 3.1. For (3 < (3 CT , we have sup u(Hp) = 0, while sup cr(Hp) = Xo(/3) > for 
j3 > j3 cr . In the latter case, Xq(/3) is a simple eigenvalue. It is a strictly increasing and 
continuous function of f3. Moreover, lim^^ Xa(fl) = and lim^oo A(/3) = oo. 

Due to the monotonicity and continuity of A = A (/3) for (3 > (3 cr , we can define the 
inverse function 

/3 = A)(A):[0,oo)^[/3 cr ,oo). 

Let C = C\(-oo,0]. 

The resolvent Rp{\) = {Hp — A/) -1 is a meromorphic operator valued function on C 
Let us introduce the operator 

A(\) = v(x)R (\) : C cxp (R d ) -+ C cxp (R d ), A e C. 

We have the following identity for the resolvent 

Rp(X) = Ro(X)(I + (3A(\))-\ A e C. (12) 

From this it is not difficult to show that l//3o(A) is the principal eigenvalue of —A(X) and 
using this, we can extend the domain of (3 to [0, oo) U (U PI C), where U is a sufficiently 
small neighborhood of zero. 

Lemma 3.2. (Asymptotic behavior) 

a) The principal eigenvalue has the following behavior as (3 I f3 CT 

V/3) = i(/3-/3 cr ) 2 (l + o(l)). 

I Per 

b) The small X asymptotics of /3q(X) is given by 

-tttt = J- -7iv / A + C(A), A-K),AeC. (13) 
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Lemma 3.3. We have 

a) For every e > 0, Rq(X) : C cxp (R 3 ) — > C(R 3 ) is a bounded operator with ||i?o(A)|| = 
C(|A| _1 ) as A -> oo, | argA| < n - e. 

b) The operator (I + PA(X))^ 1 on C cxp (IR 3 ) is meromorphic in C. For each e, A > 0, 
there is a 5 > such that it is uniformly bounded operator in A G C, | arg A | < 
7r-e,|A| > A, \0-p a \ < I 

As for the A — > asymptotics, 

Lemma 3.4. There are Ao > and 5q > such that for X G C'U{0}, |A| < Ao, — /3 cr | < 
5 ,/3 7^ j5o(A), we aai>e 

(J + /M(A))- 1 = + S(X)) + C(A, /?) 

as operators on C exp (M 3 ) ; where B is the one dimensional operator with kernel 

B ( x ^) = 7 ? \ l2 ( w ' 

J R3 f(x)^ 2 (a;)o(a; 

5 = 0(a/|A|) ; 5(0) = ; and C(X,f3) is bounded uniformly in X and (3. 



Now for / G C ex p(M 3 ), define 

gl r (z,y) = (I + P(T)A(z))- 1 f(y) 
The key to the proofs in this paper is the following lemma. 
Lemma 3.5. For 5 > 0, 

as T oo for X G T(7 2 /2 + 5), where K(X,T) is uniformly bounded as an operator on 
C exp (M 3 ). 

Proof. First consider the |A| < aT case for a sufficiently small in order for the following 
to make sense. Lemma 13.21 yields 

Po{z) = Per + PlllVi + O(Z), Z^0,zeC. (15) 

Using that 

/3(T) = /3 cr + -j= + o C±=\ , T ^ oo, (16) 



we can write for A G T(7 2 /2 + 5) that 

A, (^) - P(T) = #71^-- "J= + ° (f ) + ^V^) 
as T — 7- oo. Using the definition of 7, this can be rewritten as 

MP M|) = 



(17) 



/3 cr7l (v / 2A-7) 

where C a (l) denotes a bounded quantity depending on a. 
Combining f ll8p and Lemma [3.41 gives 

/ ( x \ _ 271 Im^( x )f( x ) dx v(y)1>(y) 



TV 

+ a(i), 



/3 cr (/ R3 ^(^(x)^) 2 a/2A - 7 



as |A| < aT and T — > 00, where O a (l) is a bounded operator. 

For |A| > aT, the first term on the right hand side of ( TT4"|) is uniformly bounded. By 
Lemma I3TB1 and /3(T) — > (3 CI , so is the left hand side, and therefore their difference too. □ 

Proof of Proposition \2.£\ Fix 5 > and note that 

u m (t,x) = ^- [ e xt (R (\)g f T (K-))(x)d\= (19) 

Zm Jr(X (j3(T))+S/T) 

If x± f e-^^7T f\ \ dX 
= / e t / . -g J \ y\ dy— 

*™ Jr(T\ (p(T))+6) JR3 2ir\x-y\ \T J T 

after a change of variables, where we used the explicit expression for the kernel of Rq(X), 
which is just Green's function for the Laplacian on the whole space. Note that by Lemma 
13.31 a) . moving the contour is permitted. 
By Lemma 13.21 and (fl~6l) , we have 

Ao(/?(T)) = -^r(/?Cn - P^) 2 + o((/3(T) - /3 cr ) 2 ) = + o (i) , (20) 

and we get TXq(/3(T)) — >■ -^W- = 7 2 /2. Therefore, for large enough T, we can take the 

Tl Per 

path T(7 2 /2 + 5) as the contour of integration. 

Let us denote the first term in ( fl4l) by g?, and let g\ = gr — g%- If u%t\ stands for 
the contour integral ( fl9l) with g^ in place of gx, then we have 

u° f){T) (t,x)= (21) 

1 «(/) /" e A ^ If e-^-^T \ 1 / \ 1 1 \ 

v{y)i){y)dydX. 



VT 2tu Jr(4+6) V2X - 7 J R3 v(w)ip(w)dw J R3 |sc - 2/ 
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It's easy to see by Taylor's formula that for y bounded, ey/T < \x\ < e 1 V / T, A G 
T(7 2 /2 + 5), there are C, a > such that for large enough T, 



e ' \/r 



e vt 



\x - y\ 



x 



< C- 



T 



Plugging this back into ( |2TT) . the first term gives the main term of (El), while the remainder 
is easily shown to satisfy (jHJ) (as v has compact support, all integrals exist). 
The remaining error term is 



u 



0(T) 



(t,x) 



2niT 



\ 1 



e 

2n\x — y\ 



(K(X,T)f)(y)dydX. 



Splitting the spatial integration as 

1 



u 



0(T) 



(t,x) 



2niT 



T(^+5) J{\y\<^f} 
we get, after making the substitution £ := Xt/T, 



h + h 



la 



1 1 

2m T J T 



Now l/\x -y\< 2/(e v / T) yields 



141 < 



2tt\x — y 



{K(X(0,T)f)(y)ey e* dyd£. 



Before we estimate h, we need to make the following observation. The contour of 
integration in all the previous formulas as well as in the expression for lb, can be bent 
towards the negative real axis. Namely, by T'(a), we mean a union of two rays emanating 
from a and that make a ±45 degree angle with the negative real axis. In all the preceeding 
formulas, the integration can be performed on either of the corresponding contours T 
or r" since the integrands are analytic and decay along the imaginary axis and decay 
expononentially in the negative real direction. 

Then the change of variables £ := Xt/T yields 



1 1 



f 

Setting x = \[Tz and y = VTu, we get 



1 



2n\x — y 



{K(X(0,T)f)(y) e y e~y dyd^. 



\h\ < 



(27T) 



y3/2 



r'(A( 7 2/ 2 +5)) J\u\> 



\z — u\ 



)T dud^ 
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from where the exponential decay of | h | as T — > oo and thus the first claim of the Theorem 
follow. 

It is not difficult to deduce ([9]) from ([7]) after noting that the fundamental solution at 
time t is the solution with the initial data pp{t, y, 5) evaluated at time t — 5. □ 
Proof of Proposition \2.3[ Let ut = Pp(T) — Po- Then 

^u T = H p{T) u T + /3vp , 
with initial condition zero. By the Duhamel formula, 

u T (t,y,x)= / pp(T)(t - s,z,x)fiv(z)p (s, y,z)dzds, 

Jo Jr 3 

which can be written as 

u T (t,y,x) = 



1 k f e T vt 



VT2m Jr(^ + s) (V2X--y)\x 



d\^(z)(3(T)v(z)po(s, y, 0)dzds + hxit, y, x) 



where h T (t,y,x) is an error term. The first term, denoted by u°(t,y,x), can be easily 
seen to equal 

1 (3(T) f f 

u%(t,y,x) = — =— — / w T>X)X (t - s)p o (s,y,0)ds, (22) 



a/T Per JO 

where 



1 r e "T-^vT 



w T ,x, x (t) = it-. / 2 r-d\. (23) 

2m Jr(^+8) (V2A -7)|x| 

We can evaluate the convolution in f )22|) in the following way. First note that (123]) is 
an inverse transform, while the transform of po(s,y,0) is e~^^\ y \/2ir\y\, and thus the 
transform of the convolution is 



T e -V2A(|rr| + |y|) 



27r( v /2AT-7)|x||?/| 

Applying the inverse formula of the Laplace transform and substituting A — > A/T, we get 
the main term in ([TO]) . 

The remainder term can be written as 

h T (t,y,x) = hft(t,y,x) + h%\t,y,x), 

where 

D 1 

lo Jm? VT 2ni Jr{^+8) (V2A~7)k 
■ip(z)v(z)/3(T)(p (s, y, z) - p (s, y, 0))dzds. 



h$\T,t,y,x)= I I -4^/ „ 7^= — d\- 
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Using the same Laplace transform trick, this can be shown to equal 
1 /3(T) 1 f f /V^V e ^ 



y/T Per 2ni J rC A +s) y/2\ - 1 \x\ J R3 \ \y -z\ \y 



dzdX, 



and this can be shown to satisfy ffTTj) the same way the main term in ([7]) followed from 
(12 ip . The last remaining term hj) containing the error term of (|9]) can easily be shown 
to satisfy ([II]). □ 
Proof of Proposition 2Jj_ Applying the Laplace transform techniques, it is not difficult to 
show (see Lemma 7.1 in [3]) that 

Z m>t (x) - 1 = ~ f ^(R m (X)((3v))(x)dX. 

Using this formula, one can prove the claim following the same steps as in the above proof 
of Proposition 12.21 □ 
The next lemma easily follows from the Feynman-Kac formula. 

Lemma 3.6. For = t < t\ < ... < t n < T, 

P/3,t(£(*i) e dx l ,... ) x(t n ) e dx n ) = — dx x ...dx n . 

Proof of Theorem \2.1\ By Lemma I3.6[ the finite-dimensional densities of the measure 
fx-^P P(t),t f° r < t% < ... < t n < 1 and x%, ...,x n G M 3 are 

Pti,...,t n ( a 'l5 • •) X n) — 

= T 3n/2 Pm (hT, 0, x 1 T^)...p m {{t n - t^T, Xn-iT 1 / 2 , x n T^) ^(^(i-M (^ 1/2 ) 

^(t),t(UJ 

Let's introduce 

p T (s, t, y, x) = p m {T{t - s),yT l l\ xT 1 ' 2 ) 

and 

T („ + „\ z fj{T),T(i-t)( xTl/2 ) 



R T (s,t,y,x)=T 3 / 2 



Pl3(T)( S ^^y^ X )z 0(T) [ ni _ a) (yT^) 



t < 1, 



pj m (s,l,y,x) j _ ^ 



Then it is not hard to show that 

pL..a( Xi ' ■■■i x n) = i? T (0,ti,0,Xi) • ... • J R T (t„_i,t n ,x n _i,x n ). 
Note that by Proposition 12.21 and Proposition 12.31 we have for i^O 

lim T Pm (Tt, 0, xT^) = ?Mf dX, (24) 



T ^ U 2vrz J rC A +5) (y/2\-y)\x 

11 



while for x, y ^ 0, 



lim T 3 / 2 p p(T) (Tt, yT 1 ' 2 , xT 1 / 2 ) = (25) 

T— »oo 



e xt-V2\Qx\+\y\) 



dX. 



(2vrf)3/ 2 AnHj TC _l +s) (V2A- 7 )Nl2/| 
By Proposition \2A\ for x 7^ 0, 

Z^t^xT 1 / 2 ) Z 7)t (x), T -»• 00. (26) 
Using this, (|2"5|) and (TJJ, it follows that for x, y =fi 0, 

1- nT/ + \ TD( 4- \ - (4. •.fs3P>y(t,l> X 1 Z )dZ 

hm i? (s,t,y,x) = R{s,t,y,x) := p {t - s,y, x) . _ -7- 

J R 3P 7 (s, l,y,z)dz 

Moreover, as follows from ( |24"|) and (j4|), it is not difficult to see that if x 7^ 0, 
lim i? T ( S , t, 0, x) = t, 0, x) := limpet - s, y, x) S^P^^ z ) d \ 

Since i?(s, t, y, x) is the transition density of the polymer under the measure Q 7 (as 
discussed in Section 1 and as shown in jlj), this implies the convergence of the finite- 
dimensional distributions and the result will follow once tightness is shown. On the other 
hand, the proof of tightness is only a slight modification of the proof of Lemma 10.5 in 
[3] (where the case of fixed (3 was treated), so we don't provide it here. □ 
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